Multiple inflation and the WMAP 'glitches' 
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Observations of anisotropics in the cosmic microwave background by the Wilkinson Microwave 
Anisotropy Probe suggest the possibility of oscillations in the primordial curvature perturbation. 
Such deviations from the usually assumed scale-free spectrum were predicted in the multiple in- 
flation model wherein 'flat direction' fields undergo rapid phase transitions due to the breaking of 
supersymmetry by the large vacuum energy driving inflation. This causes sudden changes in the 
mass of the (gravitationally coupled) inflaton and interrupts its slow roll. We calculate analytically 
the resulting modifications to the curvature perturbation and demonstrate how the oscillations arise. 
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I. INTRODUCTION 

Despite the overall success of the standard ACDM inflationary model in matching the results from WMAP, it is 
striking that the goodness-of-fit statistics for the data are rather poor In particular, Xos/'^ — 1432/1342 for the fit 
to the TT spectrum which implies a probability of only ~ 3% for the best fit model to be correct. The lack of power 
at large scales relative to the standard ACDM model prediction has motivated several alternative proposals for the 
spectrum and nature of the primordial fiuctuations ^ IE IE Qi Hi 1^ Q , however the uncertainties (and 
cosmic variance) at low multipoles are large and it has been argued that the observed low quadrupole (and octupole) 
are unlikely only at about a few per cent level ^13. .14111511 ."^ In fact attempts to judge the goodness-of-fit by eye can 
be misleading since neighbouring C;'s are correlated (being 'pseudo-Q's evaluated on a 'cut' sky); taking this into 
account it is found 1] that the excess comes mainly from sharp features or 'glitches' in the power spectrum that 
the model is unable to fit. In particular there are glitches at I ~ 120, 200 and 340 amongst the first and second 
acoustic peaks. 

The WMAP team have noted that these glitches may have been caused by beam asymmetry, gravitational lensingof 
the CMB, non-Gaussianity in the noise maps, or the method of power spectrum reconstruction from the CMB maps []|. 
However they have also considered the possibility that these glitches are due to features in the underlying primordial 
curvature perturbation spectrum [l^ . specifically in the 'multiple inflation' model [l9| where sudden changes in the 
mass of the inflaton can generate characteristic localized oscillations in the spectrum. This is physically well motivated 
since spontaneous symmetry breaking can occur during inflation for 'fiat directions' in supersymmetric theories and 
such fields are gravitationally coupled to the infiaton 19]. The resulting change in the potential of the infiaton field 



was modeled earlier as 



rally 

m 



1 + Campi tanh 



C^grad 



(1) 



which describes a standard 'chaotic infiation' potential with a step starting at 0stcp with amplitude and gradient 
determined by the parameters Campi and dgiad respectively. This was shown to result in oscillations in the primordial 
curvature perturbation spectrum by numerically solving the governing Klein-Gordon equation 0| . The WMAP team 
found that the fit to the data improves if such oscillations are allowed for, with a reduction of by 10 for the model 
parameters (/)stcp = 15.5 A/p, Campi = 9.1 x lO""* and dgrad = 1-4 x 10~^Afp 18]. 

We calculate the curvature perturbation spectrum from multiple infiation more realistically, taking into account the 
evolution of both fields as dictated by the inflationary dynamics , and using a WKB approximation |23] , rather than 
numerical integration to solve the governing equation, in order to gain insight into how such oscillations are generated. 
In a subsequent publication we perform flts to the WMAP data with the cosmological parameters unconstrained, in 
order to investigate the sensitivity of their inferred values to such deviations from a scale-free spectrum |22j . 

The possibility that the WMAP glitches are due to oscillations in the primordial spectrum has also been considered 
in the context of inflationary models invoking 'trans-Planckian' physics |2E |3 l25l l2^ ^'^^'^P features in the 
primordial spectrum may also have been generated by resonant particle production '28', '29', '30] . Several authors have 
attempted to reconstruct the curvature perturbation from the WMAP data and have noted possible features in its 
spectrum 0, |^ 112, l33L 133. Issf . It is clearly necessary to have an analytic formulation of the modiflcation to the 



^ However there is an unexpected alignment betwee n th e quadrupole and octupole llSlIld and there appear to be systematic differences 
between the north and south ecliptic hemispheres llTl . so this issue cannot be considered settled as yet. 
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usually assumed scale- free spectrum of the inflationary curvature perturbation, both in order to provide a link between 
CMB observables and the physics responsible for the glitches, and to distinguish between the various suggestions for 
such new physics. 



Supergravity (SUGRA) theories consist of a 'visible sector' and a 'hidden sector' coupled together gravitationally 
[3a.l37l |. The Standard Model particles are contained in the visible sector and supersymmetry (SUSY) breaking occurs 
in the hidden sector. Most SUGRA models of inflation have the inflaton situated in the hidden sector. This is because 
it is easier to protect the necessary flatness of the inflaton potential against radiative corrections there. The visible 
sector is usually assumed to be unimportant during inflation. 

However this might not be true if one of the visible sector fields undergoes gauge symmetry breaking. Suitable 
candidates for this are the so-called flat direction fields which generically exist in supersymmetric theories (for a 
review, see ref . |38j| ) . These are directions in field space in which the potential vanishes, in the limit of unbroken 
SUSY. During inflation in supergravity theories scalar flelds typically receive a contribution of ~ zLO{H^) to their 
mass-squared due to SUSY breaking by the large vacuum energy driving inflation USES 113. (For the inflaton this 
constitutes the notorious '77 problem' |42| : in common with most models of inflation we need to assume that some 
mechanism reduces the inflaton mass-squared by a factor of a least 20 in order to allow sufficient inflation to occur.) 
The SUSY-breaking mass term typically dominates the potential along the flat directions. If the mass-squared is 
negative at the origin in fleld space, the flat direction p will be stabilized at an intermediate scale vev of 



by higher dimensional operators (x p"/Mp jlOj. These appear in the potential after integrating out heavy degrees 
of freedom, because we are working in an effective field theory valid below some cut-off scale, which we take to be the 
reduced Planck mass: Mp = (SttGn)"^/^ ^ 2.4 x 10^^ GeV. 

We assume that the inflationary era when the observed curvature perturbation was produced was preceeded by a 
hot phase, so that all fields were initially in thermal equilibrium. The potential along the flat direction is then [ijj 



with a smooth interpolation aX p ^ T . Here we have included the 1-loop flnite temperature correction to the potential 
psLli^l with Nyi{T) being the number of helicity states with mass much less than the temperature. 

The finite temperature correction therefore creates a barrier of height 0(T'^) in between the origin and p ^ T'^ /m. 
The tunneling rate through the barrier is negligible and so p is confined at the origin until the temperature, which 
falls rapidly during inflation, reaches T ^ m and the barrier disappears. The flat direction fleld p then undergoes a 
phase transition as it evolves to its global minimum at E. This is governed by the usual equation of motion 



II. MULTIPLE INFLATION 
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which has solutions 
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(5) 



E -I- i^i exp ( 



3Ht \ ■ 3Ht fi~ o^ 8m^ 

— jsm _^(n-2)g^ 



I + K2 , {p)^^. 



Here we have taken 




(6) 



Therefore p evolves exponentially towards its minimum, with most of the growth occurring over the last ~ 2 e-folds 
[T9l |. and then performs a few strongly damped oscillations about E, as seen in Fig.^ (We assume that damping 
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e-folds 

FIG. 1: The evolution of the flat direction field p, corresponding to two different values for the change in the inflaton mass. 

by particle creation as p oscillates about its minimum is less important than the severe damping due to the large 
inflationary Hubble parameter evident in Fig.^) 

After the phase transition the vacuum energy density driving inflation is reduced by a factor of [l — (S/Mp)^] . For 
E significantly below the Planck scale this can be neglected, hence H can be taken to be sensibly constant. A more 
important effect occurs due to the gravitational coupling between p and the inflaton field 4>. 

The particular background inflationary cosmology in which the phase transition occurs is not directly relevant for 
our work. For definiteness we take it to be a 'new inflation' potential with a cubic leading term which is a phenomeno- 
logically successful model with a naturally small inflaton mass ^^^|. Then, for example, a term Kcfxp'^ p^ / Mp in the 
Kahler potential which is allowed by symmetry (near 0^0) yields a coupling term A0^/9^/2 in the inflaton potential, 
where A = Km^/Mp jlJi]. With 

V{c^,p)^Vo-C3cl)' + ^X<l>^p' + ..., (7) 

this causes the effective mass-squared rn^ = d'^V/d(p'^ of the inflaton to change from to^ — — 6c3((/)) before the phase 
transition, to — —603(0) + afterwards. This will affect the curvature perturbation spectrum which is very 
sensitive to the mass of the inflaton. 

In order to produce observable effects in the CMB or large-scale structure the phase transition must take place as 
cosmologically relevant scales are leaving the horizon. In supersymmetric theories there are many flat directions which 
can potentially undergo symmetry breaking during inflation, so it is not unlikely that a phase transition occurred in 
the ~ 8 e-folds which is sampled by observations [13 , particularly since the observations suggest that (the last period 
of) inflation may not have lasted much longer than the minimum necessary to generate an universe as big as the 
present Hubble volume. 

III. THE CURVATURE PERTURBATION 

Some care is required in calculating the curvature perturbation spectrum for multiple inflation because we might a 
priori expect two-field inflation during the phase transition, when cj) and p evolve simultaneously. This might mean 
that both curvature as well as isocurvature perturbations are created, involving a significantly more complicated 
computation ^|. However, as we show below, our model of multiple inflation is effectively single- field inflation as far 
as the generation of curvature perturbations is concerned. 

The curvature perturbation on comoving hypersurfaces TZ is given hy TZ = 6Af, where Af is the number of e-folds 
of expansion measured locally by a comoving observer passing from an initial, spatially flat, slice of space-time to a 
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final, comoving, slice |4g|. The perturbations 50fts and Spfts are defined on the spatially flat time slice and the final 
comoving slice defines IZ. The first time slice is taken during inflation, soon after horizon crossing, and the second 
time slice is taken at the end of inflation, after TZ has become constant. During the phase transition the slow-roll 
conditions in are no longer satisfied but we assume the end of inflation to be well after this epoch. Then we have 

0(p op 



and so 

^-^l^j ''^"l^J 

which is just the single-field inflation expression for V-r- This is because the vacuum energy is dominated by 0, which 
drives inflation, and p has a negligible effect on the number of e-folds of inflation. 

The curvature perturbation spectrum is usually calculated analytically using the slow-roll approximation to some 
chosen order in the slow-roll parameters |42l |. This assumes the potential is relatively smooth, which is however 
inappropriate for multiple inflation. For example, for the standard chaotic inflation potential V ~ ^rn^cj)^, the usual 
slow-roll expression gives 



1 

TZ'' — 



1/3/2 



dV/d(t) 



1 1/3/2 

(10) 



where all quantities are evaluated at horizon crossing. This gives a simple 'step' in the amplitude if the inflaton mass 
changes suddenly as in eq. Q j missing all the fine detail that we expect in the spectrum. (A discussion of the slow-roll 
approximation applied to a potential with a sudden change in its slope |4g] rather than in its curvature can be found 
in ref.^j.) 

Therefore we must resort to first principles and use the general formalism 51] for calculating the curvature pertur- 
bation spectrum. Instead of the slow-roll approximation we use the recently suggested WKB approximation [22, |53| ■ 
The metric describing scalar perturbations in a flat universe can be parameterized as 

ds^ = [(1 -I- 2/1,) drf - 2d^Bsdrjdx' - {(1 - 2Ds) Sij + 2d,djEs} dx'dx^] . (11) 

We employ the gauge-invariant quantity 

u = a(50(s'' -I- = zU, (12) 
to characterize the perturbations; see refs. fssLls^ lssI for more details. Here 



b' {B, - Ei) , (13) 



a 



are also gauge- invariant, z = a4i/H, and 



* = D,--{Bs- E'J , (14) 
a 



n = Ds + H^. (15) 

The primes indicate derivatives with respect to conformal time 

f 1 

where the last equality holds in de Sitter space. The Fourier components of u satisfy the Klein-Gordon equation of 
motion 

The spectrum depends directly on Uk and is given by 



(18) 
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The solutions of ea. (|17|l are governed by the relative sizes of fc^ and z" j z. In de Sitter space z" jz = therefore 
initially when 3> z" / z and the mode is well inside the horizon, we have the flat spacetime solution 



Uk 



2k 



(19) 



Uk = AkZ + BkZ 



(20) 



where iji is some arbitrary time at the beginning of inflation. In the opposite limit, k^ <C z" /z, we have the solution 

when Uk is well outside the horizon. The first term on the right is the growing mode solution and the second is the 
decaying mode. Thus V-ji at late times is equal to k^ |^fe|^ /27r^. 

We will use the WKB approximation to solve the mode equation, with ea. (|19|l as an initial condition. The WKB 
approximation cannot be directly applied to ea. H17() because the approximation breaks down on super-horizon scales. 
However, as noted in ref.|23|, the WKB approximation becomes relevant if we switch to the variables 



In 



The transformed mode equation is then 



where the effective frequency. 



U 



d^ 



Ha 



e^/'uk. 



Q{x)U = 0, 



Q{x) 



1 



-2x 



generally decreases with time during inflation, passing through zero at least once.^ In de Sitter space, this is 

Q{x) ^' ^ ' 



k 

In the WKB formalism |5^ we expand U as an asymptotic series 



U (x) — exp 



.71=0 



, Sn (x) > S'„+i (a;) . 



Substituting this into ea. (|23|l results in the following series of equations 

So = ^Q, 
2S'oS[ + Si; = 0, 

2-5*0 -Sn + S',"-! + X! '^j'^n-i = ^' for n > 2, 

where the primes indicate derivatives with respect to x. 
This gives the first few terms as 



Q>0 <^ 



' So 




Si 




S2 




^s. 


_ Q" 

16Q^ 



Q" I 5Q\ 
8Q3/2 32Q3/2 



(21) 
(22) 

(23) 
(24) 

(25) 

(26) 



(27) 
(28) 

(29) 



(30) 



^ Q is related to the variable gs of Refs.l53. l57l by Q = 



6 



Q<0 < 



(So = ±r(-Q)^/^dy 
Si = -iln(-Q), 

^2 



±r 






L 8(-Q)^''^ ' 32(-Q)3/2j 


Q" 


5Q^ 




16Q2 


64Q3 • 





(31) 



Using only the first two terms in the expansion and neglecting the rest gives the Ist-order WKB approximation (for 

g >o) 



A 



exp 



)l/2 



B 



■ exp 



)l/2 



(y)dy 



(32) 



where A and B are constants. At early times Q ~ e which means that the WKB solution will match ea. ll9|l if 
A = and B = l/\/2fc. 

The solutions become inaccurate close to the roots ('turning points') x» of Q(x,) — 0. To continue the solution 
through a turning point we therefore need to match the WKB solutions valid to the right and left of the turning point, 
with a solution valid in the neighbourhood of the turning point. We patch the solutions together using asymptotic 
matching. 

We consider the case of a single turning point and divide the x-axis into three regions. Region I is defined as the 
region where x < x*, Q > and the WKB approximation is valid. Region II is the neighbourhood of a;, where Q can 
be approximated by its tangent at x,. Region III is where x > x,, Q < 0, and the WKB solution is accurate. We 
assume there is an overlap between regions I and II, and also between regions II and III. 

In region II we have 



dg 

dx 



Q~—aX, X = x — X*, a. 
We rewrite the integration limits of Ui (|32|l as 



> 0. 



Then in the overlap between regions I and II, Ui becomes 

TT ^^^^^ I ^^-l/4 
Ui ~ — -^==- (— A ) ' exp 



2k 



With Q ~ — aX, the solution to ea. H23() is 

Uii ~ CAi {a^/^X 



DBi 



i (a^'^x) 



(33) 



(34) 



(35) 



(36) 



where Ai and Bi are Airy functions of the first and second kinds and G and D are constants. For A" ^ 0, Uu becomes 



-1/12 



2V^ 



(-A) 



-1/4 



{D - iC) exp { i 
2 



X 



,3/2 ^ - 

4 



+ {D + iC)c^^{-i[^{-Xf^ 



Requiring this to match Ui H35|) in the overlap region fixes the coefficients 



C = iD, D^J—a 



,-1/6 



exp 



For X we have 



-1/12 

7^ 



^-1/4 



c 



-exp ( --ai/2A3/2 ) +i^exp ( -a^'^X^' 



2 1/ 



(37) 



(38) 



(39) 
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In region III, the Ist-order WKB solution is 



F 



-Q{x)] 
G 



-Qix)] 



^exp|£[-Q(y)]i/^dy| 



Fexp + Gexp (-^a'/^X'/^ 



When Q ~ —aX, this becomes 
which matches Un (|39|) in the second overlap region if 

Substituting Um into ea. (|18|l gives finally the curvature perturbation spectrum 



3/2 



(Xf)]'^/^expi / [~Qiy)f'dy 



(40) 



(41) 



(42) 



(43) 



where X{ is the value of x at some late time well after horizon crossing when TZ is constant, and we have neglected 
the decaying mode proportional to G which is negligible at late times. 

We now consider the 2nd-order WKB approximation where the first three terms in the expansion H26(l are retained. 
At early times the WKB solution is 



Ui 



1 



2fcQi/4 
1 



2fcQi/4 



exp 



exp 



)l/2 



Q" 5Q 



'2 



)3/2 32Q3/2 



48Q3/2 J^^ 48Q3/2 



dy 
" Q" 



(44) 



where the second line follows after an integration by parts. As before we rewrite the integration limits as 



Q" 



48Q3/2 



dy = / Q'/^dy~ 
Q"'dy 



X,—6 ^11 



Q" 



3/2 

Q" 



Q'/'dy- 



3/2 

Q" 



dy 
dy, 
^dy. 



(45) 



Here e is a small parameter introduced to regularize the divergent integral 5*2 . This time, region II is the neighbourhood 
of a;* where 



-aX + PX'', (3= - 



Id^Q 



2 dx2 



(46) 



Then in the overlap of regions I and II 

rv-l/4 



2k 



(3X 



4a 



3/2 



£-1/2 + T 



, /3 



-3/2 



-1/2 



12a3/2 



(47) 



With Q ~ —aX + ea. (|23|l has the approximate solution 

Uii k(^- 

+L I 1 



^^-^ Ai 



5a J 
5a 



Bi 



ai/3x (^1 
a^/^X ( 1 - 



5a 

(3X 
5a 



For X < 0, 



t/ii 



-1/4 



/3X 
4a 



^3/2 ^ 

5ai/2 



(L - ^K) cxp ^ I ^ (-X)^'^ + (-X)'^/^ + ^ 



which matches Ui in the overlap region provided that 



exp 



4 12a3/2 



,-1/2 



When 0, we have 



4a 



-cxpf-V/'X3/2 



5ai/2 

The 2nd-order WKB solution in region III is 

M 



+ .exp(V/^X3/2-^xV^) 



Uiii 



i-Q) 



1/4 



exp 



N 



[-Q) 

M 



1/4 



i-Q) 



1/4 



exp 



5g 



'2 



8(-gf/' 32 (-Q) 



3/2 



dj/ 



Q" 



5Q 



'2 



(_Q)3/2 32 (_Q) 



3/2 



48 (-Q) 



3/2 



48 (-g) 



3/2 



dy 



dy 



N 



i-Q) 



1/4 



exp 



5Q' 



48 (-Q) 



3/2 



With Q ~ —aX + PX"^, this becomes 



x.+e 48 (-Q)^/' 



dy 



48ai/2 



^-3/2 



12a3/2 



12a3/2 



+jVexp ( --ai/2x3/2 + %- 

^ ' 3 5ai/2 48ai/2 



^-3/2 



12a3/2 



X-l/2. 



/3 



12a3/2 



-1/2 



Hence if 
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then Uiii matches Uu. Therefore in the 2nd-order approximation 



-,1/2 



ik) 



2T:(j) \aH 



3/2 



exp 



f3 



X exp 



12a3/2 
5Q' 



'i-Q) 



-1/2 



3/2 



Q" 



3/2 



dy 



(55) 



A weak dependence on the parameter e arises because of the incomplete cancellation between exp (^Y2a^^ ^^^) ^'^^ 



computational convenience 



above. This is because we have truncated the expansion of Q at 2nd-order (see ea. H46() ') for 



IV. RESULTS 



During the phase transition, the scalar potential can be written as 

V {cj,, p) = Vo- cscj,' - m^p^ + IxctPp" + -^P". 

2 Mp 

Then the change in the inflaton mass-squared after the phase transition is 



Ami = AE^ 



The equations of motion are 



n-4 ni/("-2) ^2m^Af"^''^ 



(2m2 - X^-") 

717 



n7 



dV 



(3c30 - Xp^) 0, 



p + 3i7p= -f = ( 2m^ - A<^^ - ^p"- , p. 



Using these gives 



— = ( 2i/ 2 + 6c30 - Ap^ 
z 



2App0 \ 



which is shown in Fig.[21 taking 



0.05i/^. The effective frequency of ea. (|24ll is then obtained to be 



i/2 



2 



= 7^^ I ^ - 7^' " 6C30 + Xp' + ^^^'^ ^ 



(56) 

(57) 

(58) 
(59) 

(60) 
(61) 



We calculated the WKB solutions using (j) (t) and p (t) found by numerical solution of eas. H58|) and H59|) . An example 
of a Ist-order WKB solution is shown in Fig.|31 while Fig.0] shows the power spectrum obtained by numerically 
integrating ea. (|17|l mode by mode, together with the Ist-order WKB approximation H43|l. for two different values of 
Am^. The spectrum V-r, has the form of a positive step with superimposed damped oscillations and a spectral index 

n — 1 = ''j"^^^ of n ~ 0.961. The Ist-order WKB spectra are too small by about 10% at large and small scales. They 
do exhibit oscillations but these are slightly out of phase. We can understand the spectrum as follows. 
The equation Q = can be rewritten as 



k^s{t), sit) = xl— + 



2^2 



(62) 



Therefore £ = a/s defines a new length scale in inflationary cosmology — the 'potential scale'; it isL not the horizon 
scale as is commonly stated in the literature, that controls the form of the mode functions [53|. When a mode 
is inside the potential scale it has the oscillatory behaviour of eq. (|39|l and when it is outside it has the exponential 
behaviour of ea. (|32|l . In de Sitter space we have £ = 2/3H~^ and so the two scales are nearly equal for slow- roll 
inflation. However, they behave quite differently during the phase transition, with £ undergoing oscillations. 
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FIG. 2: The evolution of z" / z during the phase transition, for two different values of the change in the inflaton mass. 
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FIG. 3: Comparison of the Ist-order WKB solution (1361 with ut calculated numerically for k — 6.7 x 10 



turning point is at a;, ~ 8.64 and Am^ 
is continuous and close to the exact solution. 



O.lml. 



Note that the solutions Ui 13211 and I7iii 14UI I diverge at 



hMpc"\ The 

X, but Uu ira 



The oscillations in the curvature perturbation spectrum occur on length scales which cross the 'potential scale' as 
it is oscillating. The maxima of the oscillations in V-jz correspond to scales which cross at the minima of the potential 
scale. These modes begin their exponential growth early and so have increased by a larger amount at X{. Similarly 
modes which cross at maxima of i give rise to minima of Vn- The amplitude of the oscillations of the spectrum 
decreases with wavenumber because the oscillations of z" / z are damped. 

Modes which cross the potential scale well before the phase transition are growing according to Ufc cx z by the time 
of the phase transition, and so Viz is unchanged on large scales. Modes which cross the potential scale well after the 
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FIG. 4: Comparison of the Ist-order WKB approximation with the numerically calculated exact spectrum, for two different 
values of the change in the inflaton mass. 



phase transition arc insensitive to the variation of z" j z at the time of the phase transition, because k'^ ^ z" j z then. 
However, they cross outside the potential scale slightly later than would be the case if there were no phase transition, 

6c3(/)] to I — V 6c3(/) — AE I . Therefore for fixed Xi the modes will 

have undergone less growth. On the other hand, the introduction of the phase transition causes an even larger decrease 
in the value of (/)f and so T'k increases on small scales. Consequently there is a step in the spectrum. The factor by 
which the phase transition causes the spectrum to increase is given using the Ist-order WKB approximation as 



because i has increased from 



3c3(0f-)" 



_ XT? ' 



X exp 




IPT 



1/2 



,1/2 



(63) 



Here the superscripts 'PT' and 'NPT' denote quantities evaluated with and without the phase transition, and we have 
used the slow-roll expression (jjf ~ ~jh^ which is applicable at late times. The WKB expression l|63|l is surprisingly 
accurate as can be seen from Fig.|3 

As shown in Fig.lHI the spectrum calculated using the 2nd-order WKB approximation has an extra dip (at k ~ 
5 X IQ-'^ h Mpc-^ in this example) and is too large by ~ 10% at small and large scales. The 2nd-order solution is 
also weakly dependent on e as seen in Fig.[7| we noted earlier that this is due to the incomplete cancellation between 



exp [j^^e 1/2) and exp [/J 



Q" 



dy 



in ea. (|55() . Overall the 2nd-order WKB approximation is in fact less 
accurate than the Ist-order approximation because the asymptotic matching is less successful. As can be seen from 
Fig.lHl the 2nd-order solution is indeed more accurate that the Ist-order one away from the turning point, but is less 
accurate close to x, where it is more divergent. Hence the 2nd-order solution is less accurate in region H, consequently 
ea. (|48|l which is matched to it there is less accurate than ea. (|36() which matches the Ist-order solution. The overlap 
between regions I and II, and between regions II and III is reduced, so we have the seemingly paradoxical result that 
the Ist-order approximation is superior to the 2nd-order one. 
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FIG. 5: Comparison of the numerical result for the factor by which Viz increases due to the phase transition with the Ist-order 
WKB expression, taking A; = 10 h Mpc~^. 




1.6 



10" 



10" 



10" 



k/ (h Mpc" 



FIG. 6: Comparison of the 2nd-order WKB approximation with the numerically calculated exact spectrum, for two different 
values of the change in the inflaton mass. 



CONCLUSIONS 



It is usually assumed that the curvature spectrum from primordial inflation is scale-free, especially when fitting 
cosmological models to data on CMB anisotropics and large-scale structure. This is indeed the case for simple models 
of the inflaton where its potential is assumed to be dominated by a (fine-tuned) single polynomial term in the region 
where the observed fluctuations are generated. However it must be emphasised that these are 'toy' models and the 



13 




0.25 



0.15 



— Error of Ist-order WKB 
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FIG. 8: The error in V^kRo (uk) for the 1st- and 2nd-order WKB solutions. 

actual physics behind inflation is likely to be more complicated [sol l60l Wi\ . In particular when one attempts to 
construct realistic models of inflation in the context of supergravity or string theory, more than one field is typically 
involved and their mutual interactions are unlikely to yield the usually assumed smooth potential . Scalar fields other 
than the inflaton can exhibit non-trivial dynamics during inflation and affect the slow-roll evolution of the inflaton. In 
particular, 'flat direction' fields, which are generic in such models, may undergo symmetry-breaking phase transitions 
and, by virtue of their gravitational coupling to the inflaton, induce sudden changes in its mass, resulting in multiple 
bursts of slow- roll inflation punctuated by such phase transitions [l9j |. 

We have calculated analytically the effect of such sudden changes in the inflaton mass on the curvature perturbation, 
using a WKB method 21] to solve the governing Klein-Gordon equation without recourse to any slow-roll approx- 
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imations. Although the analytic calculation does not yet provide an exact match to the exact numerical solution, 
it does capture the characteristic oscillations that are generated in the spectrum [23, and provides considerable 
insight into the physics behind this phenomenon. There may be scope for developing more accurate calculational 
techniques in this context 's^, '57] . It should be of course borne in mind that caution is required in using the standard 
calculational framework for the inflationary metric perturbation for the case of a time-dependent effective potential 
■ Nevertheless when the scalar fields involved are weakly coupled and the time scales involved exceed the 
Hubble time, the usual results from perturbation theory ought to be valid. 

Given that the precision WMAP data does show preliminary evidence for such oscillations we are now inves- 
tigating what the data imply for the physics of the multiple inflation model and indeed whether the estimation of 
cosmological parameters from the data is significantly affected when the primordial spectrum is not scale- free p^. 
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